Conductance length autocorrelation in quasi one— dimensional 

disordered wires 

Q\ . Klaus Frahm* and Axel Miiller-Groeling^ 

^^ , Service de Physique de VEtat condense, CEA Saclay, 91191 Gif-sur-Yvette, France. 



X5 



C/3 



cd 



a 

a 
o 
o 



X 



Abstract 



Employing techniques recently developed in the context of the Fokker-Planck 
approach to electron transport in disordered systems we calculate the conduc- 
tance length correlation function {5g{L)5g{L + AL)) for quasi Id wires. Our 
result is valid for arbitrary lengths L and AL. In the metallic limit the cor- 
relation function is given by a squared Lorentzian. In the localized regime 
it decays exponentially in both L and AL. The correlation length is pro- 
portional to L in the metallic regime and saturates at a value approximately 
'^ . given by the localization length ^ as L ^ ^. 
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I. INTRODUCTION 

In the last fifteen years two powerful approaches to calculate the electron transport prop- 
erties of quasi Id disordered wires have been developed: a method based on the nonlinear 
a model by Efetov and Larkin [|l| and a Fokker-Planck approach by Dorokhov 0, Mello, 
Pereyra and Kumar [0]. Both methods are nonperturbative in the sense that they do not 
assume that the system size is small compared with the localization length ^. Therefore 
the full range from the metallic to the localized regime can be described within these two 
frameworks. Very recently it has been shown Q that both approaches are in fact equiva- 
lent, despite their considerable technical dissimilarity. Apart from the important conceptual 
point of unifying our theoretical understanding of quasi Id wires this proof of equivalence 
enables us to choose among the available techniques the more suitable one for a given prob- 
lem. Combining the results obtained within both approaches, our knowledge about quasi 
Id transport is very advanced, some would say almost complete. We mention only a few 
important facts: Both the conductance g and its variance var(f7) have been calculated for 
arbitrary system length and for all three symmetry classes |0,§]. The probability distribu- 
tion function for the generalized eigenvalues of the transfer matrix is known exactly in the 
unitary case [^ and the corresponding n-point correlation functions are also under control 
1^. However, in spite of all this progress, there are still a few unsolved problems. One 
of the probably most important questions concerns the magnetic field correlation function 
{6g{L,B)6g{L,B + AB)). This quantity cannot be formulated within the so-called "min- 
imal" a model, i.e. the supermatrix space necessary to describe this correlation function 
is larger than in previous nonperturbative calculations. The technical difficulties to treat 
larger a models exactly have so far been insurmountable. In the Fokker-Planck approach, 
on the other hand, the question is how to incorporate the magnetic field in the theoretical 
formulation. This more conceptual problem also still awaits a solution. 

In this paper we calculate the conductance length correlation function {6g{L)6g{L + AL)) 
for all L and AL using the Fokker-Planck technique. In this case the situation is more 
favourable than for the magnetic field correlation function since the system length — unlike 
the magnetic field — is the fundamental parameter of the whole approach. Therefore it 
is possible to perform an exact calculation of a quantity that cannot be formulated within 
the framework of the minimal a model. We briefiy summarize our results: In the metallic 
limit L, AL <C ^ the length correlation function is given by a squared Lorentzian, while it 
decays exponentially as a function of both L and AL in the localized regime L, AL ^ ^. As 
L increases from the metallic to the localized limit the correlation length in AL, which is 
proportional to L for L <^ C,, saturates at a constant value equal to C,- This explicitly confirms 
the view that a strongly localized wire is composed of independently fiuctuating segments 
of size ^. In a recent paper p|, S. Feng, J.-L. Pichard and F. Zeng have argued that, in the 
localized regime, the conductance fiuctuations induced by the motion of a single impurity 
should be universal, {[In g/g']"^) ^ 1. Here, g and g' are the conductances before and after 
the impurity has moved, respectively. The main qualitative assumption necessary to arrive 
at this conclusion was precisely the decomposition into segments, which is rigorously derived 
in our paper. Therefore, the universality of the conductance fiuctuations in the localized 
regime can now be regarded as analytically established in quasi Id wires. The arguments in 
M, however, go beyond the quasi Id case and suggest, supported by numerical simulations. 



that this universahty extends also to 2d and 3d. 

The conductance length autocorrelation function is in principle accessible by experiments 
in niesoscopic nanostructures. For this, one has to fabricate a "quasi one-dimensional" 
quantum-wire with several contacts at discrete distances Lj. Using such a geometry, one 
can measure the conductance as a function of Lj (and of the magnetic field 5), where the 
disorder realization is identical for the common parts of the wire. One has to take into 
account, however, that in these multi-lead devices phase coherent scattering in the attached 
contact arms typically influences the experimental result. As usual, the statistical average 
necessary to calculate the mean conductance or the length autocorrelation function may be 
replaced by an average over a suitable range of the magnetic field. 

Our paper is organized as follows. After collecting a few important known results con- 
cerning the Fokker-Planck approach in section II, we develop the technical framework for 
the length correlation function of an arbitrary linear statistic in section III. In section IV 
we derive a (complex) expression for the special case of the conductance length correlation 
function. Sections V,VI and VII are devoted to the metallic, localized and crossover regimes, 
respectively. In section VIII we give a short summary and make some concluding remarks. 
Considerations of a purely technical nature have been deferred to appendices A to E. 

II. KNOWN RESULTS AND BASIC TERMINOLOGY 

In the Fokker-Planck approach 0J^Q-|TT[ the transmission properties of a quasi Id wire 
with A^ conducting channels and length L are characterized by A^ generalized eigenvalues 
Ai > (i = 1, . . . , A^). They parametrize the radial part []3|,[l2[| of the transfer matrix. With 



Tj = (1 + Aj) ^ the transmission coefficients, the conductance g is given by the Landauer 
formula 



N N 1 



(1) 



The multiplicative combination law for the transfer matrix yields the following Fokker- 
Planck equation [|,|[|lT| for the probability distribution p(A, t) of the wire, 

dtv{\t) = /\p{\t). (2) 

Here, A = (Ai, . . . , \n) is the collection of eigenvalues Aj, t = L/2^ the length of the wire 
measured in units of (twice) the localization length ^, and A the radial part of the Laplacian 
on the transfer matrix space, 

A = AY,dxM^ + X,)jC\)dxJ-\\) , (3) 

i 

^(A)=ni^^-^.i' • (4) 

i>j 

We have ^ = {PN + 2 — P)l, where /3 = 1,2,4 corresponds to the orthogonal, unitary, and 
symplectic symmetry class, respectively, and / is the elastic mean free path. The usual 
initial condition for (H) is given by p(A, 0) = S{X), corresponding to ideal transmission 
Tj = 1 for a wire of zero length. In the unitary case (to which we restrict our attention 



in this paper) the Fokker-Planck equation (Q) has been solved exactly by means of a 
Sutherland-transformation. Subsequently, the n-point correlation functions of the Aj (at a 
fixed value of the length L) have been calculated for arbitrary A^ and L. In the following 
we summarize some of the results in [^0 for further reference. 

The probability distribution p{X, t) can be expressed in terms of a certain "many- 
body" Green's function G{il, A; t) 

p{\,t)= [ d^ilG{ll,\;t) p{il,0) = \imG{il,\;t) . (5) 

This Green's function is given by 

G(/i, A; t) = ^ ^ det ((^(/i„ A,; t)) e^-* , (6) 

where 

P(A) = n(^.-A,) , Cjv = -E^- ' £„ = (l + 2n)2 (7) 

i>j n=Q 

and g{fi, X; t) is a "one-body" Green's function corresponding to the differential operator 

D{X) = -{4dxX{l + X)dx + l) . (8) 

We have 

dtg{fi,X;t) = -D{X)g{fi,X;t) = -D{n)g{fi,X;t) , g{fi,X;0) = 6{fi - X) . (9) 

The operator D{X) has a continuous set of eigenfunctions 

^.(A) = Pi(,,_i)(l + 2A) = F(i + ^|,l- 2|; 1; -A) (10) 

with eigenvalues /c^, i.e. D{X)'ipk{^) = ^^V^fc(A). Here, Pi(j;i._i) denotes the generalized Leg- 
endre function, which can be expressed in terms of the hypergeometric function F{a,b]C]z) 



as indicated above []T3|. The expansion of g{n, A; t) in terms of these eigenfunctions is given 
byi 

roc 2 

gifi,X;t)= J^ dA;|tanh(f)^,(^)^fc(A)e-'^* . (11) 



Using (0) and (11) the limit in (Rl) can be evaluated B. The resulting expression can be 



rewritten M in the very useful form 



p(A,t) = ^det(g„_i(A,;t)) det(/i^_i(Ai;t)) , (12) 



where (n, m = 0, 1 A^ — 1 



5 -^5 



Q„(A;t) = P„(l + 2A)e-^"* , (13) 

hM;t) = J^ dkltanh{f)LUk')M^)e-''' , (14) 



The Legendre polynomial Pn(l + 2A) in (0) is also an eigenfunction of -D(A). The cor- 
responding eignevalue is — £„. The P„(l + 2A) do not contribute in the expansion (pj]) 
because they are not normalizable in the range of integration A > 0. The advantage of the 
representation (|12D is due to the biorthogonality relation [^ 

("OO 

/ dX Qn{X]t) hm{X;t) = 5nm (16) 

which was the key to calculating the n-point correlation functions [^. 
In the following, we will also need the two properties 

/■oo 

/ dfihm{l^;t) g{fi,X;At) = h^{X;t + At) , (17) 

Jo 

/ dXg{iJ,X;At)Qn{X;t + At)=Qn{l^;t) . (18) 

Jo 



These identities can be verified using (H), (p^ ) and the fact that dtQniX;t) = D{X)Qn{X;t) 
(seeim). 



III. DENSITY DENSITY CORRELATION FUNCTION FOR THE 

EIGENVALUES 

We are finally interested in calculating the conductance length correlation function 
{Sg{L)6g{L + AL)). To this end we need to know the joint probability density function 
of finding a certain fixed set of eigenvalues fi sX t and another fixed set A at t + At. Obvi- 
ously, this density is given by 

P2(/i, t-Xt + At) = pill; t) G(/x, A; At) , (19) 

since the conditional probability density of finding a transition /i — > A in a length intervall 
At is precisely given by the propagator G{fi, A; At). 

We are now going to formulate a slightly more general problem than the one we wish to 
solve in the end. Let y4i2 = Z]jCti,2(Aj) be two arbitrary linear statistics of the eigenvalues 
Aj. Then the length correlation function of these two quantities can be written as 

/•OO /'OO 

{Ai{t) A^it + At)) = dXi d/ii ai(Ai)a2(/ii) i?i,i(/xi,t;Ai;t + At). (20) 

Jo Jo 

The function Rii{fii, t; Ai; t + At) is the probability density of finding one eigenvalue equal 
to /ii at t and one eigenvalue equal two Ai at t + At. It arises from the joint probability 
density (|T9|) upon integrating out all the remaining eigenvalues, 

-Ri,i(/^i, t; Ai; t + At) = A^^ / d^2 ■ ■ ■ dfj,N dX2 ■ ■ ■ dXN P2(/i, t; A, t + At) . (21) 

Using the fact that det(Q„_i(Aj; t)) = const. p(X)e^^^ and with the help of (||), (|T2D , and 
(pj|) we can reexpress p2 as 

P2(/i, t; A, t + At) = — det(/i^_i(/ii; t)) det(^(/i„ A,; At) det(g„_i(A,; t + At)) . (22) 



Rather than calculating all the integrations in (^T]) explicitly, we will employ a generat- 
ing functional |14] to calculate the density (0): Let Mi.2(A) be arbitrary functions of the 



eigenvalue A. Following |]T^ we define the functional 



L[ni,n2] = (j[ui{K{t)) n«2(A,(t + At))\ (23) 

= j d^'fij d^^X n('"i(/^0^2(A0) P2(/i,t;A,t + At) 

so that the correlation function ( pO]) can be expressed as 

{A^{t)A2{t + M))=d,,d,,L[l + z^ai,l + Z2a2] ^ . (24) 

Inserting (^) in (^) we find after some algebra involving transformations of the various 
determinants 

L[mi,M2] = det (Mm„[Mi,M2]) , (25) 

with the A^ x A^-matrix 

/•CXD /"OO 

Mmn[ui, U2] = d^i I dX hmifi] t) g{jj,, A; At) Qn{X; t + A) Mi(/i) M2(A) (26) 



depending linearly on ui and U2. To calculate the derivatives in ( P^ we notice that due to 
(p!6|)-([r8|) we have Mmn[l, 1] = ^nm so that (in matrix notation) 

M[l + ziai, 1 + ^2^2] = 1 + ZiM[ai, 1] + Z2M[1, 02] + ziZ2M[ai, 02] 

= 1 + X . (27) 

To expand the determinant in (|25|) we employ the relation 

det(l + X) =exp(tr ln(l + X)) = l + tr(X) + |(tr(X)2-tr(X2)) +••• . (28) 



Finally, using (|25|) , (^TQ , and ( psf ) in (|^) we can perform the derivatives with respect to the 
auxiliary variables zi and Z2 and arrive at 



{Ai{t) A2{t + At)) = tr(^M[ai, 02]] + tr(^M[ai, l]j tr(^M[l, 02]] - tr(^M[ai, 1] M[l, 02]^ 

(29) 

This result can now be straightforwardly compared with the expression (|20|) to identify the 
two-point density Ri^i{fi, t; A, t + At). To write Ri^i in a compact and appropriate form we 
define the function 

N-l 
i^Ar(A,ti;/i,t2) = XI <5m(A;ti)/l„(/i;t2) , (30) 

m=0 

in terms of which the one-point density -Ri (yU, t) reads 



i?i(/i;t) = ii'Af(/i,t;/i,t). (31) 

This quantity occurs in the calculation of simple averages of a linear statistic such as 

(A(t)) = 9,,L[l + 2;iai,l][^^^ = tr(M[ai,l]) =y'~t//iai(/i)i?i(/i;t) . (32) 



With the definitions ( pO]) and ( |3lD and taking into account the relations (|I6[)-(|18]) the com- 
parison between ( pOf ) and (p9|) results in 



i?i,i(/i, t; A; t + At) = i?i(/i; t) Ri{X] t + At) 

+Kn{X, t + At; /i, t) \g{fi, A; At) - ir^(/i, t; A, t + At) 



(33) 



To summarize the result of this section, the correlator of Ai{t) and A2{t + At) is given by 
(Mi(t) 5A2it + At)) = {Ai{t) A2it + At)) - (Ai(t)) (A2(t + At)) 

/•OO /"OO 

= / dX dfi ai{X) a2{fJ,) S{fi,t; X;t + At) , (34) 

Jo Jo 

where we have introduced the density-density correlation function for the generalized eigen- 
values of the transfer matrix, 



S{n, t; A; t + At) = Kn{X, t + At; fx, t) g{n, X; At) - Kiyif-i, t;X,t + At) 



(35) 



In the next section we derive an explicit expression for the conductance length correlation 
function starting from (Rl) and 



IV. CONDUNCTANCE LENGTH CORRELATION FUNCTION 



The correlator {5g{t)Sg{t + At)) is a special case of (p4D, obtained by choosing oi 2(A) 



{6g{t) 5g{t + At)) = / dX dfi —— — — 5(/i, t; A; t + At) 
JO Jo 1 + /i 1 + A 



(36) 



The evaluation of this expression is a little involved and makes use of a few mathematical 



identities, which we state as we go along. To begin with, inserting (|35D in (^) we have 
{5g{t)5g{t + At:)) = Ci + C2, where 



C2 



00 11 

dXd^i——--- — KN{X,t + At;fi,t)g{fi,X;At) , 
1 + A 1 +/i 

00 11 

dXdfi Kn(X, t + At; /i, t)KN(lJ^, t; A, t + At) 

1 + A 1 +yU 



(37) 
(38) 



The second contribution, C2, is the easier one to calculate. Replacing the K^ in (|38| ) by the 
r.h.s. of (3C) one realizes that the relevant remaining integral is given by 



dfi 



l + /i 



hmifJ',t)Qm'{l^,t) . 



(39) 



Using the expansion (see appendix A) 

^^ g„(l + 2A) = (-l)-e-^-* ( ^-p^ + ^ a^, e'^'Q,(\ + 2A) j 



(40) 



the abbreviation 



m 1 

aw=2(-l)'+i(l + 2/) ^ - (a^, = for / > m) , (41) 

!/ = « + ! ^ 



/■°° 1 

F^{t)= dii -^ hm{li]t) , (42) 

JO 1 +/i 



and the biorthogonahty relation ([T6|) it is not difficult to see that 
and hence 

/7V-1 \ /7V-1 



roc 1 

/ dfi —— hmifi; t) Qiifi; t) = (-l)'e-^'* F^t) + {-if aimC^'--'^^' , (43) 

Jo 1 + w 



^2 = - E (-l)'"e-^™(*+^*)F^(t) E i-lYe-^'^Fiit + At) 

\m=0 J \l=0 J 

N-1 

E i-ir^'a^ni [F^it + At) e-^'^*-^-* + F„(t) e--'"(*+^*)+-'^*) . (44) 



mJ=0 



The slightly more complicated contribution Ci can be evaluated as follows. We again replace 
Kjq by the r.h.s. of (BUI) and employ the decomposition (ffOl), so that Ci = C^ + C4 with 



^J ^^ roc 

C3= E(-ire-^'"^*^^*W c?Aci/x— — ^/i^(/x;t)^(/x,A;At) , (45) 

m=0 -^0 1 + A 1 + ^ 

^-1 roc 1 

C4 = E (-ire^''"''"^^*^^*^«™W c;/x^/i„(/x;t)gK/^;^) • (46) 

Jl^n Jo 1 + U 



To simplify ( |iB| ) we have taken advantage of (|TB|). Using the decomposition (^Dp in (^) and 
recalling (|T6|) we get 



7V-1 



C4 = E (-l)™+'a™«e^'^*-^'"(*+^*)F„(t) . (47) 



m,l=0 



Next, to evaluate C3 we insert the relations (|lTl) and (0) in ( ^5|) to get 






^3= E / ciA;i^tanh(-i)L„(A;2)(-l)-e- 



cifc2^tanh(-^)e-'=^''^*/ dfi-^^,,{fi)4j,,{fi) I dX-—■^Pk,{X) • (48) 
2 2 Jol + /i Jol + A 

For further progress we need the intgral 

8 



r a r 



as derived in appendix B for arbitrary complex a with Re(a) > 1/2 and 

(-1)^4(1 + 2m) a{N,m,k) 



(49) 



Lfnyk 



a{N,m, k) 



r(iv + i + .|) r(iv + i 



r(i-z|)r(i + z|) 



2 ^2 



r(iV-m) r(A^ + m + l) 



(50) 



which follows by a straightforward calculation from (|T5|). Furthermore we notice that the 
function Fmif) defined in (^2]) can be written by virtue of (|14D, (^9[) and (|50| ) as 



F„(t) = (-1)"4(1 + 2m) dk\ tanh(^ 



^;-,^a(Ar,m,/c)^_fc2i 



fc2 + £. 



(51) 



After all these steps {6g(t)6g{t + At)) is expressed as the sum of C2, C3, and C4 as given 
in (0), (^Hj), and (0), respectively, after taking into account the relations (^9D-(|5TD. Ob- 
viously, a simplification of the notation is called for. Let us therefore define the symbol 



]f '-) by 



[/(A:,m)]f'")= Y. dk |tanh( 



4(l + 2m)a(A^,m, A;) 



{/(fc,m)}( 



-(fc^+£m)t 



(52) 



with f{k,m) an arbitrary function of k and ?7i. With this notation the conductance length 
correlation function is given by 



{Sg{t)5g{t + At)) = Si{k, At) e 



-Smi^t 



-Sm'^t 



{k,m) 



(k,m) 



-k^At 



S'2(m, At) e 

{k,m) 



-fc^At 



(k,m) 



(53) 



with 



Siik,At) 



r(i-i|)r(i + i|)7o "^i + fxJo 



dfi 



^k{^l) 



°°^^_^(yU,A;At) 



dA;-tanh( — )^| If-^r Ire 

2 ^ 2 ^r i-^l r i + ^l 



1 + A 
I{k,k), 



m—l 
82(711, At) = J2 (^ml 
1=0 



m 1 u—1 



-1)' e'^'^* = -2^-^(1 + 2/) e 



^e,At 



u=l ^ 1=0 



and the abbreviation 



Iik,~k) 



00 ][ 

(i/i 

1 + /i 



Mi^)i^M ■ 



(54) 
(55) 

(56) 



In general, the summations and integrations in (|53|) have to be carried out numerically. The 
most difficult numerical task is the evaluation of (|56D . In appendix C we derive an expression 



9 



for I{k, k) which allows for an efficient numerical treatment. For completeness, we mention 
that the average conductance can be calculated from (152) and is given by |M] 



-i{k,m) 



{9{t)) = [ir"'^ . (57) 

In the limiting cases of a metallic wire (A^ — *> oo,t <S 1) and of a localized wire (t ^ 1) 
further analytical progress in the calculation of the conductance length correlation function 
(^) is possible. The next two sections deal with these two limits, respectively. 



V. METALLIC REGIME 

In this section we show that in the metallic limit t ^ 1 and for large channel numbers 
N ^ oo the conductance length correlation function is given by a squared Lorentzian, i.e. 

{6g{t)6g{t + At)) = ^j^^^^+0{t) . (58) 

While this result is rather simple its derivation requires considerable effort. We begin by 



observing that for A^ — > oo the coefficient a{N,m,k) (see (|50D ) tends to unity ||T3| so that 
the following decomposition rule for the symbol [. . .]\ holds: 

{-d,) [Mk) /2(m)]f '") = [Mk)f^ [f2{m)t^ . (59) 

Here, fi{k) and f2{jn) are arbitrary functions and we have introduced the independent k- 

and m- "averages" : 

[h{k)f^ = r dkktB.^h{f)e-'''{h{k)] , (60) 

oo 

[f2{m)t^ = E 2(1 + 2m) e^-* {h{m)] . (61) 

m=0 

Our derivation of (^) relies essentially on the small t-behavior of the quantities [1]^ and 
[l]j , which will shortly be seen to become important. 



|lll"" = ^(l + 5« + ^«')+0(«')+0(e-'''"*"') ■ (63) 

These expansions are derived to all orders in t in appendix D. In (^) we have also indicated 
the presence of non- analytical contributions. 

Let us first consider the case At = 0, i.e. the case of universal conductance fluctuations. 
Using (D, (In]), (H), and (^ one can show that ^i(fc,0) = \{l + k'^) and ^2(^,0) = 
—m{m — 1) = i(l — Em)- Hence we have from ([53| ) 

(fc.™) /rii(fc,m)\2 _ ir-.i(fc)r,i(m) /r-, i(fc,m)\ 2 



{5g\t)) = \[e+e^\\^ '- [[ir-^) = mtiir - ([i]r'"0 • (64) 
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With the help of the expansions (|62D, ( |63D and using (p9D we obtain 



[lf[ir^ = 4^ + ^ + ^(^') 



iWriiM 



[lt"'^ = - J dt {[iril]^ ^ ^^ ^^ 



l_±t+0(t3) 



(65) 
(66) 



(We omit a formal proof that the integration constant in ( |66[ ) is zero. A nonzero integration 
constant contradicts the fact that {5g{tY) = 0{1) in the metallic regime). Inserting (|65D 
and (|66|) in (Bl) we find the well-known result 



{Sg\t)) = - + 0{t') 



(67) 



As a side remark we mention that we have verified up to 0{t^'^) (using computer algebra 
and higher order terms in (|62|) and (|63D , see appendix D) that the corrections to (dglt)"^) 
in ([67| ) are even in t with positive coefficients. This shows explicitly that the non-analytic 
corrections indicated in (|63|) are vital for the onset of localization. 

Now we have to generalize our treatment to the case At > 0. We define a new variable 
X by At = tx and will only keep the lowest relevant order in t but all orders in x. The 
expansions of Si{k, At) and 82(171, At) in powers of At read (see appendix E) 



^i(A;,At) = ^ 



-At)" 



n=0 

00 



n\ 



-irr^^A-k' 



-At)" 
S'2(m, At) = Y^ — -1 — rn+i{emj 



n=0 



n\ 



(68) 
(69) 



where rn+i{z) is a polynomial of degree n+ 1 in the variable z. In appendix E, we calculate 



the first three coefficients a„, 6„, c„ in 



with the result 



Tn+liz) = anZ""^ + 6„2;" + CnZ 



.n-1 



+ 



(70) 



4(n+l)2 

f {2n + lf 
I2n{n + 1) 

1 
L 4 



(n>0) , 
(n> 1) 

(n = 0) 



(2n + l)(2n-l)(12n2-5) 
180(n + l)(n-l) 

5 
16 



(n > 2) 

(n = l) 
(n = 0) 



(71) 
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From (|68|)-(|70D it is clear that we will need the following type of moments 



{k 



2\n 



(k) 



[e 



rait 



{-dtT[lf^ 

{-dtnit 



Til 1 7 



n\ 1 7 , , 



0{t 



2-n\ 



(72) 
(73) 



which have been calculated using (|62D and (|63D. Combining (pSf) and (|69|) with (|7^) and 
([f3|) we obtain 



[^i(A;,At)] 



(fe) 



1 „ /At\ U fAt\ . /At\ ,^, , 
ai — ]+-bi — ]-c{ — ]+ 0{t) 



[S2{m,At)] 



(m) 



t2 V t 

1 /At 



t2 



t 



t V t 



t J \ t J 

+ c( — )+0(t) 



t 



with 



alx) 



b{x) 



cix) 



1 OO 1 

^ n=0 ''^ 



1 oo 
■^ n=0 



1 



-X) 



6(1 + x) 



^(i + i)Mi + .) 



6 



7 1^1 

60 2 ^, n 



1 7x 



4x^ 



+ 



2x^ 



1 

377 + 



10 45 15(1 + x) 15(l + 2;)2 V72 240a; 240 
With the help of the decomposition rule (pD|) we have 



7 7x 

+ — -)ln(l + x) 



:~^t){5g{t)Sgit + At)) = \[S^{k,At)] 



m 



-k^At 



(-...») u|,., 



-EmAt 



[52(m,At)]J'") 



To proceed we need the following relations 



-EmAt 






(m) 

t 



(fc) 



-k^At 



'SmAt 



(k) 



(m) 



-At 



1 1 T_ 

2(t + At) ~ 6 ^ 60 



( 

*+^* 2(t + At) ' 6 



(t + At) + 



[1] 



(m) 



+ i + ^« + A<) + 



-fe^At 



(k,m) 



- I dt 
1 



-fc^At 



(^) ^n(m) 



In M H ) ( \ At I (t + aiAt) + 

4 V t yVAt 3 30 ) 45^ ^ 



-£m.At 



{k,m) 



Jt 



1 , / At\ / 1 1 7 , 
-lnl + — -— + - + — At 
4 V t / VAt 3 30 



^(t + a2At) + ... 



(74) 
(75) 

(76) 
(77) 
(78) 



(79) 

(80) 
(81) 
(82) 



(83) 



Here, ai and 0^2 are two undetermined constants. However, adding (|8^) and (|83|) and 
expanding both sides for small At we find the constraint ai + 02 = 1. It turns out that our 
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final result does not depend on the single remaining unknown constant and we obtain from 



.o.)(smHt^m^-~j^-'^,o.[^^) . (84) 



The result (|58|) follows immediately upon integrating this last relation. In principle we 
still have to worry about an integration constant that might depend on At. However, this 
constant must vanish for At = due to (^) and possible higher order terms are at least 
0(At) = 0(t). 

VI. LOCALIZED REGIME 

In this section we derive the fact that for t = L/2C, 3> 1 (localized regime) but arbitrary 
channel number N the conductance length correlation function depends exponentially on 
both t and At. 

In the limit t ^ 1 only the term with m = and the region of small k contribute 
significantly to the expression (|5^) . This observation simplifies (^3]) enormously. Since 
[. . .]t '"^ ~ e~* we can neglect higher powers of this symbol and with 6*2(0, At) = we get 



{Sg{t) Sg{t + At)) ~ \Siik, At) e"^'"'^*] ^'''™^ ~ {5g\t)) 45i(0, At) e"^* . (85) 



Here, we have used the fact that {Sg{tY) ^ ^Ik"^ + £m\t ~ i[l]t '"" (see (|6^)) in the local- 
ized limit. The average {Sg{t)'^) can be calculated by means of a saddle-point approximation 
to give 

3/2 
(<5/(t))^a(iV,0,0)^t-^/2e-* . (86) 

For a(A^, 0, 0) we have from (^) 

^(.rnn. m + I? _^ (2iV)! (2iV - 1)! 

^ ' '^ r(Ar)r(Ar + l) 2 42^-i[(Ar-i)!Ar!]2 ' ^°' > 

The function 45*1(0, At) can in general (i.e. for arbitrary At) only be evaluated numerically. 
In the limiting cases At ^ 1 and At <^ 1 the following approximations can be derived 

r 1-f At + |(At)2 + ... (At<l) 

4Si(0,At)~J (88) 

i/(0,0)4^(At)-3/2 (At>l) 

where 1(0,0) = 14^(3)/7r^ ~ 1.705 is given by the integral ( [56| ) at A; = A; = 0. A particular 
consequence of our results is the very symmetric expression 

{6g{t)6g{t + At)) = const.{At)-^^h-^/^e-^'e-' (89) 

for t. At > 1. 
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VII. CROSSOVER REGIME 

To describe the crossover behavior of the conductance length correlation function between 
the metallic and the localized regime (in both t and At) we have to rely on a numerical 
evaluation of (0). We have already mentioned earlier that the key quantity I{k, k) (|56D is 
difficult to calculate as it stands and that in appendix C a more suitable expression (as far 
as a numerical treatment is concerned) is derived. 

In all the computations reported here we have set the channel number A^ to infinity. Let 
us define the quantity 

*<*'^*'= (Ss(t)sm v^T) ■ f^"' 

Obviously, this is the conductance length correlation function normalized to its value at 
At = and multiplied by the inverse of the squared Lorentzian ( pH] ) to compensate for the 
"trivial" metallic behavior. By definition, K{t, At) should be a constant equal to unity for 
t <^ 1 and not too large At. As At increases localization effects must show up and for large 
At ^ 1 we expect an exponential decay of K{t = const. <^ 1, At). For increasing t, on 
the other hand, we expect the compensation factor (1 + At/t)^ to become ineffective even 
for small At since K{t = const, At) should decay exponentially rather than algebraically 
for t > 1. In Fig. 1 we have plotted K{t,At) as a function of At for t = 0.1,0.2,0.5, 1.0, 
and 5.0. Clearly, our expectations are borne out. However, there is one intriguing and 
unforeseen feature: The curve K{t = 0.1, At) bears as function of At striking similarities 
to the behavior of the variance {6g{t)6g{t)) of the conductance in the unitary case. This 
observation suggests that (|58D remains true if we replace the numerical prefactor 1/15 by 
the full function {6g{At)5g{At)), i.e. 

(i9(«M9(t + A*))=<MMM^ + 0(() . (91) 

(l + f) 

We have not tried to prove this relation analytically and a numerical verification turns out 
to be difficult since we cannot go easily to t values smaller than 0.1. Therefore (^) remains 
an interesting speculation for the time being. 

Next, to characterize the distance over which two conductances are correlated, we have 
determined the point (At)i/e such that {Sg{t)Sg{t+ (At)i/e)) = {6g{t)6g{t))/e. In Fig. 2 we 
have plotted (At)i/e as a function of t. We see that the curve starts with a linear behavior 
(At)i/e ~ t for t < 1 and then saturates at (At)i/e ~ 0.5 for t :» 1. This is exactly what 
we expect from the limiting cases worked out in sections V and VI. In the metallic regime 
the decay scale is set by the system length itself, see (^^. This is also dictated by the 
universality of the conductance fiuctuations: There can be no intrinsic length scale other 
than the system length. In the localized regime, however, we have the asymptotic expression 
(^) and the situation changes. The exponential function, of course, does decay on a typical 
intrinsic scale. Since At = AL/2^ and the curve in Fig. 2 saturates at (At)i/e ~ 0.5 we see 
explicitly that the wire is divided into independent segments of size C,. 
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VIII. SUMMARY AND CONCLUSIONS 

Employing techniques that have recently been developed to calculate the transport prop- 
erties of a quasi Id wire in the Fokker-Planck approach we have derived the expression (|53D 
for the conductance length correlation function {6g{t)6g{t + At)). In the metallic regime 
{t <^ 1, At < 1, N -^ oo) we found a squared Lorentzian 

{6git)6git + At)) = ^j^^^ + Oit) , ® 

while in the localized regime (t, At ^ 1) the correlation function is dominated by exponential 
tails, 

{5g{t)5g{t + At)) = const.{At)-''/H-'"^e-^'e-' . (P 



For intermediate values of t and At we had to restrict ourselves to a numerical evaluation 
of the various sums and integrals in (|55|). Figs. 1 and 2 show the crossover from the metallic 
to the localized regime both in the typical dependence on At and in the dependence of the 
correlation width (At)i/e on t. The latter quantity is proportional to t in the metallic regime 
and saturates at (At)i/e ~ 0.5 (i.e. AL k. ^) in the localized regime. 

Eq. (|58|) generalizes the celebrated universality of the conductance fluctuations to the 
case of the length correlation function. As long as the ratio At/t is constant the absolute 
length of the system does not matter. The result also shows that any rearrangement in the 
wire (of the disorder potential, say) affects the full system since the correlation width is 
given by the system length. Fig. 2 proves that this changes drastically in the localized case. 
As the system size exceeds the localization length ^, the correlation length no longer grows 
with t but saturates at AL ^ ^. This demonstrates explicitly that the wire is decomposed 
into statistically independent segments of size ^. As a consequence, moving one impurity 
no longer changes the statistical properties of the whole wire but affects only the relevant 
segment. From this result one can deduce p that the conductance fluctuations {\hig/g'Y) 
in a quasi Id wire in the localized regime are universal and of order unity. 

Acknowledgements. We thank Jean-Louis Pichard for several helpful remarks. K.F. 
acknowledges C.W.J. Beenakker for fruitful discussions. This work was supported by the 
DFG and the European HCM program (K.F.), and a NATO fellowship through the DAAD 
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FIGURES 



< 



^ 




FIG. 1. The function K(t,At) as defined in the text versus At. The curves correspond from 
top to bottom to t = 0.1 (sohd), t = 0.2 (dotted), t = 0.5 (dashed), t = 1.0 (long dashes), and 
t = 5.0 (dash-dotted). 
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FIG. 2. The decay width (At)i/^ of the length correlation function versus the system length t. 
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APPENDIX A: PROOF OF (|^) 

From the definition (|T3p it is clear that we have to demonstrate that 

1 / 1 ™~^ \ 

P„(l + 2A) = (-1)- — — + Yl amiPiil + 2A) (Al) 



1 + A "- ' ^ ' Vl + A ,=0 / 

in order to prove (|iO|). Since Pm(l + 2A) is a polynomial of degree m in the variable 
1 + A an expansion of the form ( |XT| ) must exist. The coefficients ami can be calculated 
as follows. With the help of the orthogonality relation for the Legendre polynomials, 

l\ dxPn{x) Pm{x) = 6nm 2/(1 + 2n), wc get for the ami (x corresponds to 1 — 2(1 + A)) 

ami = i-iyil + 21) f' dx-^iPmix)-l)Piix) . (A2) 

j-1 1 — a; 

This integral vanishes for m < I. For m > I and using {Pm{x) — l)Pi{x) = {Pm{x) — 1) — 
{Pi{x) - 1) + Pmix)iPiix) - 1) we find 

am = i-iyil + 2l){bm-k) (A3) 

with 

' ,. ^"'^' - ' . (A4) 

-1 1 — X 

Applying the recurrence relation for the Legendre polynomials, (2m + l)xPm{x) = 
{m + l)Pm+i{x) + mPm-i{x) , we obtain the equation 

Tfl 
bm+1 — bm = —2 5m,0 H rT(^m ~ &m-l) (^5) 

m + 1 



which has the solution 



ra -| 

hm = -2Y.- . (A6) 



u=l^ 



Together with (^) this gives precisely the value for ami claimed in (^ 



APPENDIX B: THE INTEGRAL (g| 



In this appendix, we calculate the integral (PPj ) for arbitrary complex values of a with 
Re(a) > \. Using (|10D and formula 15.3.4 of Ref. we get 

V',(A) = F{\ - ^f, I + ^f; 1; -A) = (1 + A)(-i+^'=)/2F(| - ^f , | - zf ; 1; ^) • (Bl) 

Substituting s = (1+A)~^ and expanding the hypergeometric series the l.h.s. of (^) becomes 

h nl^ Jo''' ^' '^ -h nl^ Til-zl + a + n) "^^'^ 
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Here, (a)„ = a{a + 1) ■ ■ ■ (a + n — 1) = T{a + n)/T{a) is the Pochhainmer symbol [0]. The 
evaluation of the integral involves the /3-function |]13|. The sum in (B2) can be expressed as 



a certain hypergeometric function at the special value z = 1 and we can write for the r.h.s. 
of ^^ 

^^i/'i .•fc|^^ r/'i „•fc|^^^2 '2'2 2'2 ''gi"")-^; 



r(-i-2| + a)^ (1_^|)2 r(-l-2| + a) ^ ^_ ^^^_ 

■-T (^2 ^2' 2 ^2' 2 ^2 ~'~ '^' 

(B3) 



r(i-zf + a) ^Qn!(i-2f + a)„ r(i-zf + «) 



Using formula 15.1.20 of Ref. [0 we finally obtain the r.h.s. of Eq. (^ 



APPENDIX C: THE INTEGRAL (ID 



We derive an expression for I{k, k) defined in ( [56| ) that is suitable for a numerical eval- 
uation. Applying the transformation formula 15.3.8 of Ref. |T^ we have 



i^M = rafL, (1 + f^Y''^"''>^' F{\ - ^|, i - ^f; 1 - tk- ^) + {k^ ~k) . (Cl) 

i 1,2 + ^2/ 

We expand the hypergeometric function in ( Pl|) and replace the first ipkif^) in ( ^6]) by the 
resulting sum. For each term in the sum the /i-integration can be done by using the integral 
(^) derived in the previous appendix. After some simplification the result for I{k, k) reads: 

u. L. _ r(.fc)r(i-|(fc + fc))r(i-|(fc-fc)) 

^ (l-t(A: + fc))„(l-t(A:-fc))„ 1 

We cannot proceed along the lines of appendix B to evaluate the sum in ( |C2|) due to the 
presence of the extra factor {z + n)~^ (with z = ^ — i^). A direct numerical summation, on 
the other hand, is inefficient since the sum converges only like n""^. However, the convergence 
can be substantially improved by the expansion 



{z + ny ^^{2z + n)^ {z + n)^ {2z + n)m-i {z + ny{2z + n)r 

The coefficients Ci,{z), Am{z), Bm{z) depend on z but not on n and can be determined by a 
suitable recurrence relation. The contributions involving Cy{z) can be summed analytically 
by identifying (as in appendix B) a hypergeometric function F(.,.;.;l). The remaining 
terms with Am{z) {Bm{z)) behave as ~ 72-(i+™) (~ n^*^^+™)). Choosing m = 10 — 20 the 
corresponding sums converge rapidly. 



APPENDIX D: EXPANSION OF (||) AND (^ 

This appendix deals with the expansion (in powers of t) of 

[l]f)= r dkktanh{f)e 



'wk\ ^.-k'^t 



[1] 



(m) 



/O 

oo 



^2(l + 2m)e 





^{l+2m)^t 



(Dl) 
(D2) 



in the hmit t -C 1. Expanding tanh(^) in terms of exponential functions we obtain: 



[1] 



(fc) 



/ rffcA; l-2^(-l)"e 



?i^-(l+n)7rfc j g~fc^* 



^^ n=0 



dkke-'-^+''^^^ £ 



-D^f^ 



„2m 



ni=0 



m! 



2t 






2t 

1 



m=0 



m=0 

oo 



ra\ 



^2(m+l) 



n=0 



^ + l)2(m+l) 

L _ 2 y izir^ (1±M! f 1 _ ^ V(2^ + 2) 

)+ ^ Z^ ^1 ^2(m+l) l^ 02m+l j''^^"^ + ^^ 

(-1)^2(22 



2t 



1+E 



m=l 



imo/'o2m— 1 i 



22n 



IB 



m! 



2m 



(D3) 



In the last step we have used that C(2m') = (27r)^™|i?2m|/[2(2m)!], where the i?2m denote the 
Bernoulli numbers |]13[. With i?2 = | and B4 = — ^ we find the first three terms in (|6^). 

In order to treat the discrete sum (^^) we decouple the quadratic term in the exponent 
by a Gaussian integral: 



[1] 



(m) 



lim , — 

lim , 

1 



(is e-"'/* Y. 2(1 + 2m) e 



{l+2m){2is--q) 



m=0 



rfs e-^'/* -a 



lim r= / (is 

V^0+ 2t\/7lt J-00 



2 ysin(2s + ir;)^ 
2s 



-s^/t 



(D4) 
sin(2s + ir]j 

In the limit r] -^ 0+, the integrand has ^-function contributions at s = | n, where n 7^ is 
an integer. These lead to the non-analytic corrections ~ e~'^ " /i^^)^ which we will neglect 
in the following. The analytic corrections arise from small s-contributions in the integral. 
Expanding (2s)/sin(2s) in a power series (formula 4.3.68 of Ref. ||T^) and evaluating the 
Gaussian integrals we arrive at: 



[1] 



1 00 
^'' m=0 



it) 



-B. 



m\ 



2m 



r + 0(e 



^7rV(4t) 



(D5) 



Remarkably, the coefficients differ only by the sign (—1)™ from those in (^^). Again, the 
first three terms yield (B3). 
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APPENDIX E: CALCUATION OF THE COEFFICIENTS (|7lD 



In order to expand Si{k,At) we exploit the fact that the A-integration in (|5^) can be 
viewed as the formal application of the operator exp[—D{fj,)At] on (1 + /i)~^ (with D{fi) as 
in ®): 



/;.A .(M; At)^ - e-(^')- 1^] ^ E ^[-I^(.)]" T 



„=o - \- ' /^. 



oo /Aj-Nn 



= E^^n(/.) . (El) 

n=0 "■• 

Here, Rnifi') = [—D{fi)]"'{l + ji)^^ is a polynomial of degree n + 1 in (1 + /i)~\ i.e. we may 
write 

n 
Rnif^) = J2^n,li^ + f^P'^'^ , (E2) 

where the coefficient rn,i obeys the recurrence relation 

rn+i,i = -4/' r„,«-i + (2/ + 1)^ r,,; (E3) 



with ro,fc = (5o,fc. To calculate the coefficients a„, 6„, c„ in (|70D, we need the values of r„^„, 
r„,„_i and r„,„_2. From (p3D we find 



n ^l2 



-4)"n! 



n-l 

rn,n-i = (-4r-M^-l)!'E(l + 2^)' ' (E4) 

71 — 2 1/ 



>z»2 



rn,n-2 = (-4)"-^ (n - 2)!^ ^(1 + 2z/)2 ^(1 + 2i 

(E5) 



On the other hand, we may insert (^if ) and (^21) in (^) and perform the /i-integration 
using (^91), giving 



Si{k, At) = Y. ^-J-T.rn,i jj— ^,11 {(-2 + ^-^1) il + '^ + ^l)) ■ (E6) 

The product is a polynomial of degree (/ + 1) in k"^, of which we need the first three terms: 

ri((i+z.-.|)(i+.+^i)) =-Ln(p+(i+2.r 



4' 
with 



1 - ((A;2)'+i + A, {ky + Bi {ky-' + ■■■) (E7) 
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Ai = Y.{l + 2uf , Bi= Yl (l + 2z/)2l + 2z>)2 . (E8) 

u=0 Q<v<v<l 



Together, (Pq ) and (|E7| ) imply 



^^ ^ 4'^+l(n+ 1)!2 ^"'^ "^^ "^ 4n^!2 '""'"-I ' (^^) 

(-1)"+^ (-1)" (-1)"-^ 



for the desired coefficients. Combining (pif), ( p8|) and (p9|), one indeed obtains after some 



(computer-) algebra the result (|7lD . It remains to show that 5*2 (m, Ai) is given by (^) with 
the same coefficients (0). Expanding ( ^5]) we find: 



^^(m, At) = Y. ^—T^ (-2) a(m) (ElO) 

with 

CnM = E-E(i + 20'"''" • (Ell) 

u=l '^ 1=0 

For small n the C„(m) can be calculated with the help of standard mathematical formulas. 
For our purposes, however, we need the large-m behavior for all n. Therefore we consider 
the generating function 



g[m,x) 



m -| !/— 1 oo ™l+2n 



^-5:sinh[(l + 20x] = $:-— — C„(m) . (E12) 

u=l '^ 1=0 n=0 y^ "•" ^^^ 



Using (twice) the finite geometric series we can derive the inhomogeneous differential equa- 
tion 

sinh.{2x) g{m, x) + (cosh(2x) — 1) ^^.(^(m, x) = cosh[(l + 2'm)x] — coshx, (E13) 

from which we get the recurrence relation 

Cn{m) = -^ (e"^' -l-J2bnk Ckim)] (E14) 



with 



_ n + k + 2 (2n + 3\ 2(n-fe+i) .p-, .^ 



Eq. (Pi)14|) indeed gives after some algebra 

21 



(-2)C„(m) = a„£:^+^ + bnsl + CnS''^' + ■■■ (E16) 

with an, bn, Cn as in (^). 

We conclude this appendix with a short remark: Equations ( ]68| ) and (^) suggest that 
the functions 5*1 (A;, At) and 5*2 (m, At) are connected by the analytic replacement — /c^ — »• 
^m = (1 + 2r?7,)^ and At -^ —At. We have verified this for the lowest (n < 2) polynomials 
rn+i{z) and — as far as the highest coefficients On, 6„, c„ appearing in (|70D are concerned — 
for all n. This is sufficient for the purpose of the present paper. A general proof, however, 
has yet to be given. 
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